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We study the finite-size scaling of moments of the magnetization in the low-temperature phase of the two- 
dimensional Ising model. 



1. INTRODUCTION 

Recently, finite-size scaling (FSS) at first-order 
phase transitions has attracted considerable in- 
terest. Its practical importance in connection 
with Monte Carlo studies has become increasingly 
clear, for instance through the controversy about 
the order of the phase transition in quenched 
QCD. At the same time, rigorous developments 
have been accomplished by Borgs, Kotecky and 
Miracle-Sole (|]|]. These results were obtained 
for a class of models that includes the Ising model 
at large (3 and g-state Potts models with large 
q, and were proposed to be more general. The 
FSS model suggested by the results is quite pow- 
erful. It provides, for instance, an estimator of 
the transition coupling that should have expo- 
nentially small finite-size corrections. Numerical 
studies of Potts models support this prediction 
lU^], but have also found less convincing agree- 
ment for some other quantities Qg. To find out 
more about the applicability of the FSS model 
we have performed a numerical study of the two- 
dimensional Ising model Q . 

For this model, the results of ref. [Q show that 
the partition function in a large periodic box L d , 
at large /3, is well approximated by 



Z(h,L) 



,L d U(h) 



(1) 



where h = [3H, H being the external magnetic 
field. The functions f±(h) are related to the free- 
energy density f(h) — limL^oo L~ d In Z(h, L) by 

* Presented by A. Irback 



/+(_>(*) > f- (+) (h) = f(h) if h > 0(< 0), 
and /+(0) = /-(0). Each term in eq. (|l|) repre- 
sents fluctuations about one of the two coexisting 
phases, whereas the omitted remainder accounts 
for mixed configurations and is exponentially sup- 
pressed in L. This decomposition is the starting 
point for the finite-size analysis of moments of 
the order parameter. One proceeds by expanding 
the functions f±(h) about h = and obtains, at 
fixed hL d , expansions of the moments in powers 
of L~ d . The first-order finite-size corrections, to 
combinations of moments of up to fourth order, 
have been investigated numerically for the q-statc 
Potts model with q — 7, 8 and 10 in refs. As 
already alluded to, the agreement was not perfect, 
indicating that the asymptotic L regime had not 
been reached. First-order corrections have also 
been analysed numerically for the Ising model 
0, within the framework of the double-Gaussian 
model |7]||. In this earlier phenomenological FSS 
model, the two peaks in the probability distribu- 
tion for the order parameter were approximated 
by Gaussians. This corresponds to a quadratic 
approximation for f± in the model discussed here, 
as can be seen by inverse Laplace transformation. 
As a result, the double-Gaussian model produces 
the same first-order expressions, except for prob- 
lems with the relative normalization of the peaks. 
For the Ising model the normalization problem is 
absent thanks to the global spin-flip symmetry. 

We shall extend these considerations in two di- 
rections: we shall study the size-dependence be- 
yond the first-order level and we shall include in 
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the analysis moments of higher order. 

2. HIGHER ORDERS 

We consider the behaviour of the moments at 
the phase transition. We calculate FSS expres- 
sions from eq. (Q), assuming that the contribu- 
tion from the omitted tunnelling term is expo- 
nentially suppressed. For dcfinitcncss, wc restrict 
ourselves to the Ising case. The power-like vol- 
ume dependence arises from fluctuations about 
the two stable phases and is obtained by expand- 
ing the functions f± about h = 0. We write 

oc 
fc=Q 

and, by symmetry, the corresponding coefficients 
for /_ then are (— 1) dfc. The coefficient a\ = tuq 
is the spontaneous magnetization and /3a,2 = X 
the susceptibility. In order to obtain the moment 
(m k ) at h = it is sufficient to consider terms up 
to fcth order in the expansions of f±. Inserting 
these expansions into eq. (Q), it is easily found, 
for even fc(> 0), that 



k 



L- d e k2 {^)a 2 a k ' 2 



L- 2d x 



a 3 a k 3 



2„fc-4 
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(3) 



where dij = 1 if i > j and otherwise. The coeffi- 
cient in front of L~^ d is a function of a± , . . . , a J+ i , 
and it is zero if j > k. Odd moments, of course, 
vanish becuase of the spin-flip symmetry. 

For moments of moderate order k it is straight- 
forward to extend eq. (^) so as to include all 
the power-like corrections. The validity of the 
resulting expressions will however depend on k. 
In general we expect that the (approximate) va- 
lidity gets restricted to larger L with increasing 
k. Note that higher coefficients in eq. @ can be 
large not only because of the combinatorial fac- 
tors, but also because they involve higher a^. The 
series eq. is asymptotic only 0. 

Having obtained the prediction, we next dis- 
cuss how to apply and test it. For the two- 
dimensional Ising model, the spontaneous mag- 
netization <zi is exactly known, and fairly precise 



estimates of higher a^'s exist in the literature. 
Thus, one has a parameter- free prediction for low- 
order moments. In general, the a^'s are, however, 
unknown. We may then try to use fits of numer- 
ical data to eq. (||) to extract them. By check- 
ing for consistency values obtained from different 
moments, a test of the FSS model is obtained. 
In practice, such a program is not easy to follow, 
and we are therefore led to look for other possi- 
bilities. An obvious one is to try to use single- 
phase cumulants. We define the measurement of 
the single-phase cumulant by introducing a cut 
between the two peaks in the probability distri- 
bution for the order parameter. In the remaining 
part, say m > m*, where |m*| < mo, we ex- 
pect contributions from the disfavoured phase, as 
defined by the decomposition of eq. ([I]), as well 
as from mixed phases, to be exponentially sup- 
pressed in L. Up to exponentially small errors, 
we therefore expect 



a i. 



. = L { - k -^ d (m k )* 



(4) 



where the * indicates expectation value over m > 
m* . In our calculations we take m* = 0, but in 
principle one could choose any |m*j < mo- Note 
also that the applicability of eq. (H) depends on 
k. The measurement of (m fc )* rapidly becomes 
more difficult with increasing order k, since it is 
0{L~( k ~ lS>d ) and is the result of cancellations be- 
tween 0(1) numbers. 

3. NUMERICAL RESULTS 

We now would like to test the finite-size ex- 
pression for moments, eq. (|^), and also the fea- 
sibility of a numerical determination of the co- 
efficients dfc through eq. (Q). To this end we 
have performed numerical simulations of the two- 
dimensional Ising model. These were done at 
u/u c = 0.9, where u = exp(— 4/3) and u c = 
3 — 2^/2. At this temperature, the exactly known 
values of the correlation length and the sponta- 
neous magnetisation are £ « 4.8 and mo ~ 0.84, 
respectively. In addition, several of the coeffi- 
cients ak are known to high precision from power- 
series expansions ||. 

Using the Swendsen-Wang cluster algorithm, 
we simulated L 2 lattices with 16 < L < 60. For 
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each size, 4 million iterations were done. Mea- 
surements were separated by 10 iterations and 
errors were estimated through a jackknife proce- 
dure. 

The first three (normalized) single-phase cu- 
mulants Ck — L( fe_1 ) d (m k )* are shown in fig. 
They were obtained by folding the probability 
distribution onto m > 0. From the above dis- 
cussion, we expect that Cfc approaches at an 
exponential rate for large L. While our data at 
large L are consistent with such a behaviour, it is 
clear that the smaller L studied do not belong to 
the asymptotic region. Unfortunately, the data 
points left are insufficient for a study of the form 
of the asymptotic L dependence. On the other 
hand, the close agreement found at large L with 
the values in ref. Q shows that good estimates 
of the infinite volume limit can be obtained from 
the range of L studied. For instance, we note that 
the first and second cumulants are only a few per 
cent off their infinite volume values if L/£ > 5. It 
would therefore be possible to determine the L~ d 
term of the moments to a similar accuracy. 

In fig. |^ we compare our results for the first 
three even moments with the prediction of eq. (||) , 
using values of the a^'s from ref. ||. For each 
moment all non- vanishing terms in the expansion 
have been included. The curve labelled j is the 
result obtained by including the terms of order 
up to L~i d . At large L, the prediction is in very 
good agreement with the data, and the agree- 
ment is improved by adding more terms. The 
extent of this asymptotic L region already shows, 
for these three moments, a dependence on the or- 
der of the moment, as discussed above. Larger 
L is required for the sixth moment. For the sec- 
ond and fourth moment the difference between 
measured and predicted finite-size correction is 
less than 10% down to the smallest size studied, 
L/£ w 3.3. The largest discrepancy for the second 
moment, at L/£ » 3.3, is 6% or 3.4c 



4. SUMMARY 

We have studied the finite-size dependence 
of moments of the magnetization in the two- 
dimensional Ising model. The FSS model by 
Borgs, Kotecky and Miracle-Sole gives a very 



Figure 1. Normalized single-phase cumulants Ck 
against L~ d . The horizontal lines are the coeffi- 
cients at, as given by ref. [9], 
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good description of the large-L results. Data 
for the first three even moments showed, in fact, 
no significant deviation for L/£ > 6, after tak- 
ing higher-order corrections into account. The 
description remains approximately valid at the 
smaller L/£ > 3.3 studied for the second and 
fourth moment, but not for the sixth. Exist- 
ing analytical results allowed us to perform these 
tests without using any free parameters. In gen- 
eral the parameters of the FSS model are un- 
known and we therefore tried to measure them 
numerically through single-phase cumulants. We 
found that the first-order correction to the mo- 
ments could be obtained to an accuracy of a few 
per cent from lattices with L/£ > 5. 
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Figure 2. Moments (m k ) against L d . The curves 
are explained in the text. 



